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Abstract. The low energy spectrum of finite size metallic single-walled carbon nanotubes (SWNTs) is 
determined. Starting from a tight binding model for the p z electrons, we derive the low energy Hamiltonian 
containing all relevant scattering processes resulting from the Coulomb interaction, including the short 
ranged contributions becoming relevant for small diameter tubes. In combination with the substructure 
of the underlying honeycomb lattice the short ranged processes lead to various exchange effects. Using 
bosonization the spectrum is determined. We find that the ground state is formed by a spin 1 triplet, 
if 4n + 2 electrons occupy the SWNT and the branch mismatch is smaller than the exchange splitting. 
Additionally, we calculate the excitation spectra for the different charge states and find the lifting of 
spin-charge separation as well as the formation of a quasi-continuum at higher excitation energies. 

PACS. 73.63.Fg Nanotubes - 71.10. Pm Fermions in reduced dimensions - 71.70. Gm Exchange interactions 



1 Introduction 

Single walled carbon nanotubes (SWNTs) have remark- 
able mechanical and electronic properties. They represent, 
at low enough energies, an almost ideal realization of an 
one-dimensional (ID) electronic system with an additional 
orbital degree of freedom. Due to this ID character the 
proper inclusion of the Coulomb interaction between the 
, electrons in a SWNT is mandatory. For metallic SWNTs 
of infinite length the theoretical works [1,2] showed that 
I correlations between the electrons can be described within 
■ the Luttinger liquid picture. The accompanying occur- 
rence of power-laws for various transport properties could 
indeed be observed experimentally [3,4]. The effects of the 
forward scattering part of the electron-electron interac- 
tions in finite-size SWNTs were treated by Kane et al. in 
[5] within the bosonization framework. There the discrete 
energy spectrum of the collective spin and charge excita- 
tions was derived. The bosonization method has recently 
been used also to determine the transport properties of 
finite size metallic SWNT quantum dots [6]. 

So far the effect of non-forward scattering parts of the 
Coulomb interaction has only been discussed for SWNTs 
of infinite length by renormalization group techniques [1, 
2]. In [1] deviations from conventional Luttinger Liquid 
behaviour have been found only for very small tempera- 
tures T < 0.1 mK provided that the interaction is long 
ranged. The work of Odintsov et al. [2] additionally took 
into account the situation at half filling where the for- 
mation of a Mott insulating state was predicted. In the 
works treating electron-electron interactions in finite size 
SWNTs within the bosonization formalism, the effect of 



non-forward scattering parts of the Coulomb interaction 
has been neglected. This approximation, which we will call 
"standard" theory in the following, is valid if moderate to 
large diameter tubes (> 1.5 nm) are considered as in [5,6], 
or if finite size effects can be neglected since the relevant 
energies exceed the level spacing of the SWNT as in the 
experiments [3,4]. Recent experiments [7-9] however have 
found exchange effects in the ground state spectra of small 
diameter tubes which can not be explained using the "stan- 
dard" bosonization theory for interacting SWNTs. Oreg et 
al. [10] have presented a mean-field Hamiltonian for the 
low energy spectrum of SWNTs including an exchange 
term favouring the spin alignment of electrons in differ- 
ent bands. The values for the exchange energies observed 
in the experiments agree well with the mean-field predic- 
tions. However, the question of a singlet-triplet ground 
state is beyond the mean field approach. Moreover, in 
contrast to the bosonization procedure it can not predict 
the strong energy renormalization of the charged collective 
electron excitations. 

In this article we go beyond the mean-field approach. 
We derive a low-energy Hamiltonian for finite size metal- 
lic SWNTs, which includes all relevant short-ranged in- 
teraction processes. This allows us to identify the micro- 
scopic mechanisms that lead to the various exchange ef- 
fects. Using bosonization we determine the spectrum and 
eigenstates of the SWNT Hamiltonian essentially exactly 
away from half-filling. An interesting situation arises near 
half-filling since there additional processes become rele- 
vant which can not be considered as small compared to 
the dominating forward scattering terms. Unfortunately 
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we have not found a reliable way of diagonalizing the 
Hamiltonian in that situation so far. 

Concerning the ground state properties, we find under 
the condition of degenerate or almost degenerate bands, 
a spin 1 triplet as ground state if 4n + 2 electrons occupy 
the nanotube. This is insofar remarkable as a fundamen- 
tal theorem worked out by Lieb and Mattis [11] states 
for any single-band Hubbard model in ID with nearest- 
neighbour hopping that the ground state can only have 
spin or 1/2. However at the end of their article they 
explicitly pose the question whether ground states with 
higher spin could be realized in ID systems with orbital 
degeneracy, which in the case of SWNTs is present due 
to the substructure of the underlying honeycomb lattice. 
Our findings answer this question with yes, hence proofing 
that the theorem by Lieb and Mattis can not be general- 
ized to multi-band systems. Moreover it is interesting to 
notice that all of the processes favouring higher spin states 
in SWNTs involve non-forward scattering with respect to 
the orbital degree of freedom. On the experimental side 
an exchange splitting in the low energy spectrum of the 
4n + 2 charge state has indeed been observed [7-9]. How- 
ever, all the experiments demonstrating exchange splitting 
were carried out for SWNTs with a large band mismatch 
such that the ground states are supposed to be spin sin- 
glets. Especially Moriyama et al. have proven that this is 
the case in their experiments [9] by carrying out magnetic 
field measurements. Thus the threefold degenerate spin 1 
ground state has not been observed yet, since its occur- 
rence requires a band mismatch that is small compared 
to the exchange energy. Additionally to the ground state 
properties of metallic SWNTs we have also determined 
the excitation spectra. We find that the huge degeneracies 
as obtained by only retaining the forward scattering pro- 
cesses are partly lifted and the spectrum becomes more 
and more continuous when going to higher energies. Fi- 
nally this leads to a lifting of the spin charge separation 
predicted by the "standard" theory. 

The outline of this article is the following. We start in 
Section 2.1 by briefly reviewing the low energy physics of 
noninteracting electrons in finite size metallic SWNTs. In- 
cluding the Coulomb interaction we derive the effectively 
one-dimensional Hamiltonian for the low energy regime 
in Section 2.2. The subsequent examination of the effec- 
tive ID interaction potential in Section 2.2.1 allows us to 
sort out the irrelevant interaction processes. The remain- 
ing processes are either of density-density or non-density- 
density form. The former ones we diagonalize together 
with the kinetic part of the Hamiltonian by bosoniza- 
tion in 3.1. Using the obtained eigenstates as basis we 
calculate the corresponding matrix elements for the non- 
density-density part of the interaction with the help of 
the bosonization identity of the electron operators, Sec- 
tion 3.2. In Section 4 we calculate the ground state and 
excitation spectra by diagonalizing the Hamiltonian in- 
cluding the non-density-density processes in a truncated 
basis and discuss the results. 

For the hurried reader we propose to skip the more 
technical sections and, after reading Section 2, to go di- 




Fig. 1. The graphene lattice with its sublattice structure. 

rectly to Sections 4.1 and 4.2 where the results of this work 
are presented for the low energy and low to intermediate 
energy regions, respectively. 



2 Low energy Hamiltonian of metallic finite 
size SWNTs 

As shown in [12], correlation effects in metallic SWNTs are 
universal at low energies, i. e. they do not depend on the 
chirality of the considered tube. Therefore we can, without 
loss of generality, focus on armchair nanotubes from now 
on. 

In this section we will give a short summary of the elec- 
tronic structure of noninteracting finite size armchair nan- 
otubes in the low energy regime following our earlier work 
[6] . On this basis we are going to include the Coulomb in- 
teraction between the electrons, leading to an effective ID 
Hamiltonian. The subsequent examination of the effective 
ID interaction potential will determine all the relevant 
scattering processes, which are either of density-density 
or non-density-density form. 



2.1 The noninteracting System 

Before considering the effect of the electron-electron inter- 
actions, let us recall the most important facts about non- 
interacting electrons in finite size armchair SWNTs. Since 
SWNTs can be considered as graphene sheets rolled up to 
cylinders, the bandstructure of SWNTs is easily derived 
from the one of the p z electrons in the graphene honey- 
comb lattice, see e.g. [13]. Two carbon atoms p = ± occupy 
the unit cell of graphene, cf. Fig. 1, leading to a valence 
and a conduction band touching at the two Fermi points 
F = ±K e x . Quantization around the circumference of a 
SWNT restricts the set of allowed wave vectors, leading to 
the formation of subbands. For metallic SWNTs, only the 
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gapless subbands with linear dispersion, touching at the 
Fermi points, are relevant at low enough energies. Impos- 
ing open boundary conditions along the tube length L, the 
eigenfunctions of the noninteracting Hamiltonian Hq are 
standing waves (p rK (r) where the occurrence of the branch 
or pseudo spin index r = ± is a consequence of the double 
occupancy of the graphene unit cell. Furthermore k mea- 
sures the wave number relative to the Fermi wave number 
Kq and is subject to the quantization condition 



K=-(n K + A), n K eZ, \A\ < 1/2. 



(1) 



The parameter A has to be introduced if there is no integer 
n with Kq — irn/L, where L is the tube length, and is 
responsible for a possible energy mismatch ea between 
the r = + and r = — electrons. In general A depends 
also on the type of the considered SWNT [14]. Explicitly, 
¥V/e(f) can be decomposed into contributions from the two 
sublattices p = ±, 



^ = 7f E fvr (e lKX ^ P K (r) - e~ iK x Vp - Ko {r)) 



The coefficients f pr are given by 



fpr 



1/V2, p = 
-r/y/2, v = 



(2) 



(3) 



and the functions ip v p describe fast oscillating Bloch waves 
on sublattice p at the Fermipoint F, 



1 



R 



X(r-R- 



(4) 



where Nl is the total number of lattice sites and x(r 



R 



p, see Fig. 1. 

In Fig. 2 we show the linear dispersion relation for the 
standing waves ip rK . The slopes of the two branches are 
given by rhvF, with the Fermi velocity vf ~ 8.1 • 10 5 m/s. 
Including the spin degree of freedom, the Hamiltonian of 
the noninteracting system Hq therefore reads 



H = Hv F ^2 r E KC ™ 



(5) 



= + 




Fig. 2. The energy spectrum of a noninteracting metallic 
SWNT with the two branches r = ±. The level spacing is 
denoted £o and ea is the energy mismatch between r = + and 
r = — . 



By defining the slowly varying ID electron operators, 



we obtain with (2), 

«? CT (r) = VZ^sgn(F)VviMz) E ./Ww(r). (7) 

rF p 

2.2 The interaction Hamiltonian 

In this section we examine the interaction part of the 
Hamiltonian. After introducing an effective ID interaction 
potential, we discuss which of the scattering processes are 
of importance. We start with the general expression for 
the Coulomb interaction, 



is the p z orbital localized on site R of sublattice V — 



d 3 r 



d 3 r'^l(r)^,(r')U{r-r')^{r')K(r), 



where U(r — r') is the Coulomb potential. For the actual 
calculations we model U(r — r') by the so called Ohno 
potential which interpolates between Uo, the interaction 
energy between two p z electrons in the same orbital for 



r = r' and 



for large values of \r — r'\. Mea- 



4-7reoe|r— r '\ 

suring distances in units of A and energy in eV, it is given 
by [15] 



where c r(TK annihilates an electron in the state \ip rK ) \a). U(r-r') = U /\fl + (U e \r-r'\ /14.397) 2 eV 
Thus the level spacing of the noninteracting system is » 
given by 



(8) 



so 



L 



(6) 



In the next section we are going to express the interac- 
tion part of the Hamiltonian in terms of the 3D electron 
operators, which expressed in terms of the wave functions 
(p rK {r) read 

^ = E E Vr K (r)c rrTK =: "M^. 



A reasonable choice is Uq = 15 eV [16]. The dielectric 
constant is given by e w 1.4 — 2.4 [1]. Reexpressing the 
3D electron operators & a (r) in terms of the ID operators 
ifjrFa(x), cf. equation (7), and integrating over the coor- 
dinates perpendicular to the tube axis, we obtain, 



V = 



^E E sgn^^Fg^) fdxf dx'U [r][F] (x,x') 

**' {[r],[F]} 

X i J l 1 F 1 a( X )^t 2 F 2 a'( X ')^r 3 F 3 a'(x')^r i F i a(x), (9) 
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where S{[ r ] [f]} denotes the sum over all quadruples [r] = 
{r\,r 2 ,r 3 ,ri) and [F] = (Fi, F2, F 3 , F4). Under the as- 
sumption, justified by the localized character of the p z 
orbitals, that the sublattice wave functions ip pF (r) and 
P-pF{r) do not overlap, i.e., ip pF (r)ip- pF (r) = 0, the ef- 
fective ID Coulomb potential U[ r ][ F ](x,x') is given by, 

U[ r ][F]{x,x') = L 2 j d 2 r ± J d 2 r' ± ^2f pri f p , r2 f p , r J pri 

p,p' 

x <p; Fl (r)^;, F2 (r')^ p , F3 (r')^ pF ,(r)U(r - r'). (10) 

Using relation (3) for the coefficients f pr and performing 
the sum over p, p' , we can separate f/[ r ][F] into a part de- 
scribing the interaction between electrons living on the 
same (intra) and on different (inter) sublattices, 




,-- 4 





Fig. 3. The relevant scattering processes. For- 
ward/back/Umklapp scattering are denoted by f ± /b/u. 
The index I represents one of the three degrees of freedom 
r,F,a (branch, Fermi point and spin, respectively). 



U[ r ][ F ]{x,x') = - 



U$ ra (x,x')(l + ri r 2 r 3 r 4 ) 



+ U^ t ] er (x,x')(r 2 r 3 + r 1 r i ) , (11) 



where 



U™\ ra,inter {x,x')=L 2 



d 2 r ± d 2 r'. 



x l PpF 1 ( r ) l P*± P F 2 ( r ') l P±pFs{r')(PpFAr)U(r - r 1 ). (12) 

Note that the 3D extention of the considered SWNT enters 
the effective ID interaction potential via equation (12). In 
Appendix A we show how we actually determine the values 

p .1 . ■• i T jintra/ inter / r\ 

for the potentials (7 j F ] (x,x ). 



2.2.1 The relevant scattering processes 

Not all of the terms in (9) contribute to the interaction 
because the corresponding potential C7[r][F] vanishes or has 
a very small amplitude. In order to pick out the relevant 
terms, it is convenient to introduce the notion of forward 
(/)-, back (b)- and Umklapp (u)- scattering with respect 
to an arbitrary index quadruple [I] associated to the elec- 
tron operators in (9). Denoting the scattering type by 
Si we write [I]g I= f± for [I,±I,±I,I]. Furthermore we 
use [I]sj=b for [I,— 1,1,— I] and [7]s /=u is equivalent to 
[1,1, —I, —I], cf. Fig. 3. Keeping only the relevant terms, 
the interaction part of the Hamiltonian acquires the form, 



y= E E E v * 

S r =f,b,uS F =r,bS (7 =f 



(13) 



where 



V Sr s F s„ ■= ^ J J dxdx'U [r][F] {x,x') 

{Hs.[f]s F ,Ms„} 

X ' t Pl 1 F 1 a( X )^t 2 F 2 a'( X ')^r 3 F 3 a'(x')A i F i a(x), (14) 

as we are going to demonstrate in the following. 



Scattering of r We start with the possible scattering 
events related to the pseudo spin r. From (11) we can 
immediately read off that the interaction potential Um r^j 
does not vanish only if r 2 r 3 — r\r 4 . Thus we find the 
following cases for the relevant scattering types, 

i) n = r 4 , r 2 = r 3 and ii) r x = —r±, r 2 = -r 3 . 

Relation i) summarizes all the forward scattering pro- 
cesses with respect to r and the associated interaction 
potential is, 



U [r] f [F] (x,x') 



U$ ra (x,x') 



Tjinterr l\ 
U [F] \ X > X , 



:U+ F] (x,x'). (15) 



Case ii) includes all S r — b and S r = u processes and here 
the interaction potential is proportional to the difference 
between U intra and U mter , 



U [r]h/u[F] (x,x) 



U$ ra (x,x') 



U$ er (x,x') 



:Uf F] {x,x'). (16) 



Scattering of F The determination of the essential scat- 
tering processes with respect to F can be achieved by ex- 
ploiting the approximate conservation of quasi momen- 
tum. Looking at expression (4) for the wave functions 
PpF^), we find that the interaction potential {7[ r ][ F ], cf. 
(10), contains phase factors of the form e - l ( F i- F *) R x x 
e ~i{F 2 ~F 3 )R :c _ Although we are considering a finite sys- 
tem, therefore not having perfect translational symmetry, 
after the integration along the tube axis in (9) , only terms 
without fast oscillations survive 1 . The corresponding con- 
dition is given by 



F l - F A + F 2 - F 3 = 0, 



(17) 



1 For a perfectly translational invariant ID system it holds 
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that means only the Sf = f and Sf = b terms survive. 
We have explicitly checked that due to the discrete na- 
ture of the SWNT lattice also the Sf — u processes have 
very small amplitudes and can be neglected. Note that 
condition (17) leads to sgn(FiF 2 F 3 F 4 ) = 1 in (9). 

Scattering of a It is clear that only S a = f processes 
are allowed, since the Coulomb interaction is spin inde- 
pendent. 

Altogether this proofs equation (13). 



Processes conserving or not conserving the fermionic 
configuration From the discussion in Section 2.1 we al- 
ready know that we have to distinguish between electrons 
with different spin a and pseudo spin r. In the follow- 
ing we will denote the number of electrons of a certain 
species by N r „ and we will refer to the quantity TV = 
(N + ^,N + i,N_^,N_i) as fermionic configuration. Not all 
of the scattering processes in (13) conserve N. In more de- 
tail, for terms with (S r ,S a ) — (u,f + ), (S r ,So-) = (b, f~) 
and (S r ,S a ) = (u,f~) N is not a good quantum num- 
ber as can be easily verified by using equation (14). In 
general, only processes described by the N conserving 
terms are sensitive to the total number of electrons in the 
dot. As example we mention the charging energy contri- 
bution proportional to N 2 , N c := ^2 rcr N ra arising from 
the (S r , Sf, S a ) = (/,/,/) processes appearing later on. 
On the other hand for the N non conserving terms, only 
the vicinity of the Fermi surface is of relevance. 



Processes only relevant near half-filling Away from half- 
filling we find that terms with 

nFi + r 2 F 2 - r 3 F 3 - r 4 F 4 ^ 0, (18) 

i.e., the Umklapp scattering terms with respect to the 
product rF 2 can be neglected in (9) . For the N non con- 
serving terms fulfilling (18) this is a consequence of the ap- 
proximate conservation of quasi momentum, arising from 
the slow oscillations of the ID electron operators in (14) 
which near the Fermi surface are given by the exponential 

e -i[(r 1 FiN Tiai -r i F i N ri< , i )x+(-r2F 2 N r2< , 2 -r;- i F 3 N T3a:i )x'] ^fter 



/ dx' dxU(x-x')e ikx e ik ' x ' = 
Jo Jo 

F dx' f L X dyU(y)e iky e i(k+k ' )x ' = U k F dx'e l(k+k ' )x ' , 

Jo J -x 1 Jo 

where Uk = * dyU (y)e %ky does not depend on x' because 
we have assumed translational invariance. So it is clear that 
the double integral vanishes unless k + k! w 0. 

2 There are simple rules for determining the scattering type 
S r F if S r and Sf are known. Defining a product by S rF = 
S r S F = S F S r it holds, Sf + = S; S 2 = /+; f~u = b; f~b = u 
and ub — f~ . 





u+d 


ufd 




e = 1.4 


0.22A 


0.14A 


0.22A 


e = 2.4 


0.28A 


0.22A 


0.28A 



Table 1. The dependence of the coupling constants u + , uf 
and ttjf on the tube diameter d and on the dielectric constant 
e. 



performing the integrations in (14) this leads approxi- 
mately to (18). The N conserving terms obeying (18), 
Vf-bf and Vbf-f+, which describe not only processes near 
the Fermi level, add a term proportional to the number of 
electrons above half-filling to the Hamiltonian, therefore 
just giving rise to a shift of the chemical potential. 



2.2.2 Long ranged vs. short ranged interactions 

Except of U[ r ] f [F] f = ^[f] ' a ^ relevant interaction po- 
tentials £q>][F] can effectively be treated as local inter- 
actions: In the case of U^ b this is due to the appear- 
ance of phase factors e l2F ^ R!C ~ R ^ in (12), arising from 
the Bloch waves ip p F(r), cf. equation (4), oscillating much 
faster than the electron operators tp r <jF(x). The potentials 
, being proportional to the difference of the inter- and 
intra-lattice interaction potentials, are in general short 
ranged, since U^ ra {x, x') and U^ er (x, x') only have con- 
siderably differing values for \x — x'\ < ao with the next 
neighbour distance ao = 0.142 nm of the carbon atoms 
in the SWNT lattice [1]. Summarizing, only the processes 
with (S v ,Sf) = (/, /) are long ranged. All other terms 
can effectively be written as local interactions. I.e. for 
(S r , Sf) 7^ (/, /) we can use the approximation 

\ U [r]s r [F]s F { x , x ') ~ Lu Sr s F S(x- x'), (19) 
where we have introduced the coupling parameters 

u SrS F ;= 1/(2£2) J J dxdx i U[r]sr[F]sF ^ x >). (20) 

Using the approximation (19) we obtain from (14) in the 
case (S r , Sf) ^ (f,f) the following expression for the non 
forward scattering interaction terms, 

VSrSpStr ~ LU Sr S F ^ 

f L 

X / dxipl F (x)^ F ,(x)lpr 3 F 3 a'(x)A 4 F 4 a(x). (21) 
JO 

In the following we use the abbreviations u + := Uff, and 
u s F : ~ u bs F — u uSf- F° r details about the calculation, 
see Appendix A. We find that in general the coupling con- 
stants u + and Ug F scale inversely with the total number of 
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lattice sites, i.e., like 1/Ld, where d is the tube diameter. 
From a physical point of view this is due to an increas- 
ing attenuation of the wave functions for a growing system 
size. Therefore the probability of processes mediated by lo- 
cal interactions is proportional to 1/Ld. Because the level 
spacing of the noninteracting system £o scales like l/L, 
cf. (6), the products u + d/eo and uf F d/eo are constants. 
The corresponding numerical values for different dielectric 
constants e, cf. equation (8), are given in table 1. 



2.2.3 Density-density vs. non-density-density processes 

The interaction processes can be divided into density- 
density terms, easily diagonalizable by bosonization [17], 
and non-density-density terms respectively. It is clear that 
the forward scattering interaction Vf / / is of density-density 
form, 

^/// = ^EEE / / dxdx'U+ ]f (x,x') 

rr' FF> aa 1 J J 

X PrFa{x)p r 'F'a'(x') 7 (22) 

where the densities p r Fa {x) are given by 

PrFa(x) = 4>t Fa (x)lprFa(x)- 

But since we treat the short ranged interactions as local, 
also Vf+ b f+ , 

Vf+ bf + = 

Lu+ ^2 / dxipl FrT (x)ipl_ Fa (x)ip rFrT (x)ip r ^ F(7 (x) 

r*F J( > 

= -Lu + ^2j dxp rFa (x)p r - FrJ (x), (23) 

r*F J ° 

and similarly V b f+ / b f+ , 

V bf+/bf+ = 

~ Lu f+/b^2 / dxp rFrT (x)p- r±Fa (x), (24) 



3 Expressing the SWNT Hamiltonian in the 
eigenbasis of H + V pp 

Away from half-filling, the interaction is dominated by 
Vf f f . Together with Hq it yields the "standard" theory for 
interacting electrons in SWNTs [1,2,5]. Using bosoniza- 
tion we will in the next step diagonalize Hq + V pp . Sub- 
sequently we will examine the effect of V npp by calculat- 
ing the matrix elements of V npp between the eigenstates 
of Hq + V pp . The diagonalization of V npp in a truncated 
eigenbasis of Hq + V pp , discussed in Section 4 then yields 
to a good approximation the correct eigenstates and the 
spectrum of the total Hamiltonian H . 



are density-density interactions. In total the density-density 
part of the interaction is given by 

v pp = v f f i + v f+ b f+ + V b f + f + +V bbf+ . (25) 

The remaining terms are not of density-density form and 
are collected in the operator V npp . Including only the con- 
tributions relevant away from half-filling, we obtain, 

V npp - V f+ bf -+V b f+ f - + V bb f - + V uf - j + V ubf . (26) 

Near half-filling additionally the processes 

Vf-bf,V bf - f and V uf+f -, (27) 

satisfying condition (18), contribute to V npp . Overall, the 
SWNT Hamiltonian acquires the form, 

H = Hq + Vpp + V npp . 



3.1 Diagonalizing Hq + V f 



pp 



By introducing operators creating/annihilating bosonic ex- 
citations we can easily diagonalize Hq + V pp as we show 
in this section. It turns out that the Fourier coefficients 
of the density operators p r <jF{x) are essentially of bosonic 
nature. In detail, we get by Fourier expansion, 



PrFAx) = ^-Y, elSSD(F)qX P^ 



(28) 



where q — j^n q , n q <G Z. Then the operators b (TqT defined 
by, 



b<jq r ■= —=Praq r , Qr ■= T ■ q, q>0 (29) 

fulfill the canonical commutation relations \b a q-,b\, q ,} = 
Sa'aSqq' as shown e.g. in [17]. For completeness we give 
the explicit expression for b, 



aq r 5 



Joq r 



1 t 



/II, 



The bosonized expression for Hq is well known [6], 



Hq = 



q>0 



(30) 

Here the first term describes collective particle-hole exci- 
tations, whereas the second term is due to Pauli's prin- 
ciple and represents the energy cost for the shell filling. 
The third term accounts for a possible energy mismatch 
between the bands r = ±, given by 

e A = sgn(A)e mm(2\A\,\2\A\-l\). 

The operators M ra count the number of electrons N r<T 
in branch (ra). Bosonization of V pp can be achieved by 
inserting the Fourier expansion (28) into expressions (22), 
(23) and (24), thereby making use of definition (29). We 
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obtain, 



and 



v pp = v fff + v f+bf + + Vbf+/bf+ = 

(bavq + b^) 



q>0 



+ 



- u + ^2 {i> ar . q b ar . q + h.c.) 

rcr 

~ U f^2 ( b <?r-q b v-r-q + h.C.) 

TO 

~ U b E ( b °r-qbi- r . q + h.c)) \ 
ro ) 



(31) 



where the coefficients W q determine the interaction strength 
of Vf f f and are given by 

W q = J ' dx J dx'U^ F ^(x,x')cos(qx)cos(qx'). 

The last line of (31) describes the contribution of V pp to 
the system energy depending on the number of electrons 
in the single branches (re). Here E c = Wq is the SWNT 
charging energy, W c = ^2 ra N ra counts the total num- 
ber of electrons. Spin alignment of electrons with different 
branch index r is favoured by the term proportional to 
J/2 := Wj +uf. Finally the term coupling with u + coun- 
teracts the energy cost for the shell filling in equation (30) . 

Since the bosonic operators appear quadratically in 
(30) and (31) we can diagonalize Hq + V pp by introducing 
new bosonic operators a,js q and a^ Sq via the Bogoliubov 
transformation [18] given below by equation (33). We ob- 
tain 



H + V pp = ^2^2e jSq a\ Sg a jSq + \e c M 2 c 



j6 <?>0 
J 



-^■N'-ra + (e - u + ) M ra + re A 



(32) 



The first term describes the bosonic excitations of the sys- 
tem, created/annihilated by the operators Oj gq / ajs q . The 
four channels jS = c+, c— , s+, s— are associated to total 
(+) and relative (— ) (with respect to the index r) spin 
(s) and charge (c) excitations. The decoupling of the four 
modes j5, the so called spin-charge separation, will be 
partly broken by V npp - The excitation energies 6js q and 
the relation between the new bosonic operators aj$ q and 
the old operators b aqr are determined by the Bogoliubov 
transformation. In detail, we find with Sa q := £on q , 



£c+q = £o q \Jl + 8W q /e , 
£ s /c-q = eo 9 (l - u t l £ o) 



£ s + q = £o<?(l + u b /e ). 

The energies of the c+ channel are largely enhanced com- 
pared to the other excitations because of the dominating 
Vf f f contribution. For small q the ratio g q :— £Q q /£ c + q 
is approximately 0.2, whereas for large q it tends to 1 [6]. 
Small corrections due to the coupling constants and 
u + have been neglected. For the transformation from the 
old bosonic operators b aqT to the new ones cij$ q we find 

W = Yl A ™ { B J^ q + D jSg a\ Sq ) , q>0 (33) 



is 



where 




jS = c+, c— , s+, s— 

m = + t,+i,- t,- r 



(34) 

The transformation coefficients Bj$ q and Djs q in the case 
of the three modes jS = c— , s+, s— are given by 



BjSq = 1 and D jSq = 
and for jS = c+ we obtain 

VTq 



(35) 



(36) 



with g q — j^-- Small corrections to (35) and (36) re- 
sulting from the terms Vf+ & f+ and VJ, /+ /& /+ have been 
neglected. 

The physical meaning of the fermionic contributions in 
(32), depending on the number counting operators, have 
already been discussed subsequently to equations (30) and 
(31) respectively. 

An eigenbasis of Hq + V pp is formed by the states 



\N,m) :=n 

jSq 



"jSq 



|TV,0) 



(37) 



where |TV, 0) has no bosonic excitation. Remember that 
the fermionic configuration TV = (iV_|, TV_^, iV+f, TV + jJ 
defines the number of electrons in each of the branches 
(m). In the following we will use the states from (37) as 
basis to examine the effect of V npp . For this purpose we 
evaluate in the next section the corresponding matrix ele- 
ments using the bosonization identity for the ID electron 
operators. 



3.2 The matrix elements (Nm\V npp \ N'm'} 

Generally, due to V npp , the quantities TV and m are not 
conserved. Especially, the terms with S r = b, u in (26) 
mix states with different TV. However, denoting 



N s :=^sgn(a)TV r(T , 
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and 



N~ :=]Tsgn(r)iV ro 



and the bosonic part reads 

M {l] {m,m',x) = (mlc-^i^e-^iWe-^^e-^W 



e i<t>\ 3 (x) e i<pl 3 0)g l </4 4 O)gi0l 4 O) \ m > 



(44) 

we find that (JV C) N s , N~ mod 4,N~ mod 4) is conserved, In order to improve readability we have replaced the in- 
i.e, states differing in those quantities do not mix, such dices rFa by a single index L As we demonstrate in Ap- 
that the corresponding matrix elements of V npp are zero. pendix B , the explicit evaluation yields 
Note that in contrast to the real spin S z = ^N s , the 
pseudo spin S z = \N~ is not conserved in general. 

We already know that all the processes Vs r s F S a con- 
tained in V np p are effectively local interactions, i.e., of the 
form (21). Hence, in order to calculate the corresponding 
matrix elements (Nm \Vs r s F S a | N'm') we first derive an 
expression for 



M [r][F][a] {N,N',x) = 

$N , N> +E [r] w T NN , [r] [(T ] Q NN '[ r } [F](x), (45) 



where E\, 



The Klein 



M [r][F][a] (N,m,N',m',x) 



Nm 



i ; t 1 aF 1 ( X )i J t 2 a'F 2 (x)i>r 3 *> F 3 (z)VV 4 <tF 4 (x) 



N'm 1 



j • ^ri(T * &r 2 (j ! ^r 3 cr' ^r 4 <r- 

factors in (39) lead to the sign factor Tjvjv' [r] [<r] which is 
either +1 or —1 and QNN'[r][F](x) yields a phase depend- 
ing on TV". Explicit expressions can be found in Appendix 
B, equations (65) to (68). 

For the bosonic part of M^pv^-dN , ra, N' , ra', x) the 
(38) calculation in Appendix B leads to 



For this purpose we express the operators %p r<J p(x) in terms 
of the bosonic operators b ar . q and btr-q, 1 > 0, using the 
bosonization identity [17], 



M [r][F][a]{m,m' ,x) = C [r][F][a] (x) 

x A SrF (x) Y[ F(X>r][F][*] (*)> m jSg ,m' jSg ). (46) 



ipr*F(x) = ■q r<T K rcrF (x)e l ' p -^ {x, e 



i<j>l aF ( x ) e i<t>rvF(x) 



(39) 



Here the function F(X,m,m') stems from the evaluation 

e -A* a t e Aa m ,\ where 



The operator r\ ra is the so called Klein factor, which an- of matrix elements of the form ( 
nihilates an electron in the (ra) branch and thereby takes \ 

care of the right sign as required from the fermionic anti- the bosonic excitations |m)_ are created by the operators 



a\ i.e., |m) = (a^) m j\pm\ |0) . For the explicit form of 
F(X, to, to'), and the coefficients A ^upi , , (x) , see Appendix 
Vrcr \N,m) = (—1) 1=1 1 \N — e rCT ,m) , (40) g function C f [r][P , ][o-]( a; ) 1S conveniently considered in 

combination with QNN'[r][F]( x ), namely the product 



commutation relations, in detail 

,(r<r)-l 



where we use the convention I = + T,+ J.,— |, — 1= 
1,2,3,4. K rcrF (x) yields a phase factor depending on the 
number of electrons in (ra), 



K raF (x) 



■II. 



sgn(F)(r-I\f ra + A)x 



Finally, we have the boson fields i(j) r aF(x), 

< * In 

q>a 



(41) 



(42) 



K N[r\{F][a]{x) ~ Q NN> [r] [F] ( x )C[ r ] [p] [ CT ] (x) 

can be reexpressed as 

K N [r][F][a](x) = QN[r]lF](x)Cs r S F S a {x), 

where 



(47) 



In Appendix B we are going to demonstrate with the help exp 
of the bosonization identity (39), that the matrix elements 
from equation (38) factorize into a fermionic and a bosonic 
part, 



QN[r][F]( x ) — 
- 4 
3 



*7 \J2 1 = 1 S S n ( r J F 3) N r 3 ^ +^2 S S n ( r 3 F j) 



3=3 



^,4 



M[r][F][a](N,m,N',m',x) = 

M [r][F]{<r](N N', x)M [r][F][a] (m, ra', x), 
where the fermionic part is given by 

M[i](N , N' , x) = 
(N\ Kl(x) V lKl(x)rilK l3 (x)rn a K u (x)r, u \N>) (43) 



Here J2i=i a i denotes the sum ai + a 2 
Cs r s F s lr (x) we obtain 



«3 



04. For 



C f+bf -(x) = -Cf- bf (x) = -C bf - f +(x) = 

l/4sin 2 (|a;), (48) 

-C uf - f+ (x) = 4 sin 2 (^.t) 



ubf+ 



(x) 



(49) 



and Cs r s F s„{x) = 1 for the remaining processes of V npp . 
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0.4 0.6 
x/L 



Fig. 4. As rF =u(x) as a function of x for a (6,6)-SWNT. Note 
the large magnitude of As rF =u{x) compared to As rF ^ u (x) = 1 
for the processes only relevant away from half-filling! 



The function As rF (x) is differing from 1 only for terms 
with S r F — u, i.e., for the terms fulfilling condition (18) 
and which hence are relevant only near half-filling. The 
reason for this is that only for the SVf = u terms the co- 
efficients ^w[F] w W re l a ted to the charged c+mode are 
not vanishing. Hence A u (x) depends strongly on the en- 
ergy dispersion of the c+ mode and therefore on the for- 
ward scattering part of the interaction, in detail 



A u (x) = exp 



g>0 



£ c+q 



sin 2 (qx) 



Since for the repulsive Coulomb interaction eo q /e c+q < 1 
holds, we find As rF=u (x) > 1. In Fig. 4 we show As rF=u (x) 
for a (6,6) SWNT. It is the large magnitude of A u {x), 
that poses problems for properly treating the situation at 
half-filling. Moreover we can expect that even for large 
diameter tubes, interaction processes with S r F = u can 
not be neglected near half-filling. Altogether, we get with 
equations (21), (45) and (46) for the single contributions 

to V np p, 



(Nm\V Sr s F sJN'm') = 

^ ^iV,JV'+B [r]o .„/7jvs r S . 

{Hs r ,[F]s p ,Hs„} 



X — u 
AL 



N=4n: 



N=4n+\: 



N=4n+3: 



HH :::::: 
r-- 



N=(«, «,«,«) 



N=(»+l, n,n,n) 
+permutations 



fi=(n+\,n+l,n+\,ri) 
+permutations 



Fig. 5. The lowest lying eigenstates of Hq + V pp without 
bosonic excitations for the charge states N c = 4n, N c = 4n + 1 
and N c = 4n + 3. On the right side the fermionic configurations 
are given. We use the convention N = (iV+f, iV+i, iV_f, iV_^). 



exemplarily the proper calculation for (iVm | Vf + bf- \ Nm'} . 
We here give the regularized result for m = m' , since it 
is of special importance for the discussion of the ground 
state spectra away from half-filling, 

(Nm \V f + b f - | Nm) = u + ^ 

r 

i _i_ 



AL 



E 



{[r] /+ ,[F] i) ,[<r] / _} 



1 n^^wl^M^)'^'^)^'^) 1 ) • ( 51 ) 



4 The SWNT spectrum 

In Section 3.1 we wave diagonalized Hq + V pp and in Sec- 
tion 3.2 we have determined the matrix elements of V npp 
in the eigenbasis of Hq + V pp , Away from half-filling the 
magnitude of V npp is only small compared to Hq + V pp 
and therefore we can easily analyze the effect of the non- 
density-density interaction V npp on the SWNT spectrum 
by representing the total Hamiltonian Hq + V pp + V npp in 



^^AT[r][F](^)^s rF (^)n F ^[rHF]H( x )' TO i'59' m ^9)- a truncated eigenbasis of Hq + V Pl 

(50) 



4.1 The low energy spectrum away from half-filling 

We start with the examination of the ground and low en- 
ergy states. As basis we use the lowest lying eigenstates of 



of electrons N c . 



The evaluation of (50) causes no problems except for the 
N conserving processes with (S r , Sf, S a ) = (f + ,b, f~), 
if~,b,f), i b J~J + ), since then we find K N[r][F][a] ~ 

1/4 sin 2 (fx), cf. equations (47) and (48), causing the in- H o + V pp without bosonic excitations with a given number 

tegral in (50) to diverge for J2 3 s q m jS q — m 'jSq - 1> sucn 
that the evaluation of the corresponding matrix elements 
needs special care in this case. The origin of this divergence 
lies in the fact, that if no bosonic excitations are present, 
the N conserving processes depend on the total number of 
electrons in the single branches (compare to the fermionic 
contributions to Hq + V pp in (32)). Since the bosoniza- 
tion approach requires the assumption of an infinitely deep 
Fermi sea [17] this leads, without the correct regulariza- 
tion, necessarily to divergencies. In Appendix C we show 



4.1.1 N c = An, N c = 4n+ 1, N c = 4n+ 3 

First we consider the charge states iV c = An, N c = An + 1 
and iV c = An + 3. In that case the lowest lying eigen- 
states of Hq + V pp , shown in Fig. 5, which are of the form 
\N, 0) and therefore uniquely characterized by N, do not 
mix via V npp . That means that the only correction from 
Vnpp to Hq + Vpp stems from the N conserving process 
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Vf+ b f- . For states without bosonic excitations, equation 
(51) yields, because of F(X, 0,0) = 1, 

(AT, \V npp \ N, 0) = (AT, \V f+ b f - \N,0) = 

u + ^2mm(N rh N rl ). (52) 



Hence here V npp yields an energy penalty for occupying 
the same branch r. This effect has already been found in 
the meanfield theory of Oreg et al. [10]. The parameter 
5U there corresponds to our constant u + . The energies 
of the lowest lying states for N c = An, N c = An + 1 and 
N c = An + 3 only depend on N. In detail we find with 
(32) and (52), 



E N = \e c M 2 c +U+J2 min(AV T , N H ) 

r 



-—N-rcr + (SO - U + ) Mra + TE A 



(53) 



From (53) it follows that for the states depicted in Fig. 
5 the interaction dependent part of En is the same for 
all fermionic configurations N corresponding to a given 
charge state iV c . Hence the interaction leads merely to a 
common shift of the lowest lying energy levels for fixed 

Nr. 



4.1.2 N c = An + 2 

Of special interest is the ground state structure of the 
N c = An + 2 charge state, since here the lowest lying six 
eigenstates of Hq + V pp without bosonic excitations, de- 
noted |-/V, 0) with N = (n + l,n+l,n,n)+ permutations, 
mix via V npp , leading to a S = 1 triplet state and to three 
nondegenerate states with spin 0. For e^kO (the mean- 
ing of i=a will become clear in the following) the triplet 
is the ground state. In the following we are going to denote 
\(n + l,n+l,n, n),0) by |T4, -), \(n + l,n, n,n + 1), 0) by 
|t, I) and analogously for the remaining four states. Ignor- 
ing interactions, the six considered states are degenerate 
for ea = 0. As we can conclude from (32) the degeneracy 
of the six considered states is already lifted if including 
only the density-density interaction V pp , since then the 
energy of the spin 1 states |T, t) and ||, j) is lowered by 



J/2 := 



A , A 
U f +U b 



(54) 



relatively to the other ground states. Let us now con- 
sider the effects of V npp . The diagonal matrix elements 
(TV, \V npp \ N, 0) are again determined by equation (52), 
leading to a relative energy penalty for the states |||, — ) 
and | — ,TI)- Mixing occurs between the states |T,I) and 
||, t) via V b f+ ^-and V bb f - and between |||, — ) and | — , ][) 
via V u f- f- and V ub f- . With equation (50) we find 



(U \V npp \ l,T)=- 



J 



(U,-\Vn P p\-,U) 



In total, the SWNT Hamiltonian H = H a + V pp + V npp 
restricted to the basis spanned by the six states |T,T), 
11,4), IT, I), II, T), IU,-) and |-,TI) is represented by 
the matrix, 

H = E Q A n+2 + 

(-i o 





2 U 



V o 



u+ - e A 
j 

2 



(55) 



£A J 



where E 0An +2 = \E c N^ + {2n 2 +2n+l) (s - u+)-J(n 2 + 
n) + 2u + n. Diagonalizing the matrix in (55), we find that 
its eigenstates are given by the spin 1 triplet 

|T, T), |T, T), i/y/2 (|T, 4) + 14, t», 

the spin singlet 

i/V2(ira>-u,T» 

and the two states 

(c 1/2 |T4,->±|- T4», 



+ 1 

where the coefficients Ci/ 2 are given by 



Cl/2 = 



^E A + (J/2fT£A 
J/2 



Relatively to ^0,471+2, the corresponding eigenenergies are 
— J/2 for the triplet states, J/2 for the singlet state and 
u + ± \J e\ + (J/2) 2 for the remaining two states. Thus 
under the condition 



J/2 > 



{J /2? 



i.e., for a small band mismatch ea ^ J/2 the ground state 
is degenerate and formed by the spin 1 triplet, otherwise 



by 



( c 2 |T4, — ) + | — , T4))- The ground state spectra 



for the two cases ea — and ea S> J/2 are shown in 
Fig. 6 for a (6,6) armchair SWNT (corresponding to a 
diameter of 0.8 nm). Assuming a dielectric constant of 
£ = 1.4 [1], the calculation of the coupling parameters 
according to Appendix A yields values of J = 2(uf + 

u^) = 0.09eo and u + 0.03eo which agree well with the 
experiments [7,9], where nanotubes with ea 3> J/2 were 
considered. To our knowledge, experiments in the regime 
£a ^ J/2 demonstrating exchange effects have not been 
carried out so far, such that a validation of our predictions 
for this case, namely the existence of the ground state spin 
1 triplet and the mixing of the states |T4, — ) and | — , T4) ls 
still missing. The latter effect could be of relevance for the 
understanding of the so called singlet-triplet Kondo effect 
[19] in SWNTs. 
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Fig. 6. Low energy spectrum of a (6, 6) SWNT for the charge 
state N c = 4n + 2. (a) In the case ea = the ground state is 
formed by the spin 1 triplet (— > ©) and the states |TI, — } and 
— ,11) mix (— >0 states), (b) For ea S> J/2 the ground state 
is given by the spin state \H, — ). 

The spin singlet state l/\/2 (| j, J,) — \l,1)) is indicated by ®. 
The coupling parameters are J = O.OSteo and u + ~ 0.03eo. 



It should be stressed that all exchange effects, lead- 
ing amongst others to the spin 1 triplet as ground state, 
result from S r ^ / interaction processes. In the work of 
Mattis and Lieb [11] however, there is no such additional 
pseudo spin degree of freedom. Hence we suspect that this 
is the reason why their theorem can not be applied in our 
situation. 



4.2 Excitation spectra away from half-filling 



1 - 

- 




+v +v 

npp pp 

+v,„ 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

Number of Eigenvalue 



Fig. 7. The excitation spectrum for a (6,6) SWNT occupied by 
N c — 4n electrons. In grey we show the spectrum as obtained 
by diagonalizing the Hamiltonian of the standard theory H s t = 
Ho + Vf f f and in black for the full Hamiltonian H = Ho + V pp + 
Vnpp. A band mismatch ea = is assumed. The energy of the 
lowest c+ excitation is 4.3en- All other interaction parameters 
are as in Fig. 6. Arrows indicate eigenenergies of the "standard" 
Hamiltonian H at = Ho +Vf / / involving excitations of the c+ 
mode. 
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Fig. 8. The excitation spectrum for a (6,6) SWNT obtained 
by diagonalizing the full Hamiltonian H = Ho + V pp + Vn PP for 
N c = 4m and ea = 0.3eo- The spectrum becomes quasicontin- 
uous at relatively small energies. Shown are the lowest 10000 
eigenenergies. 



Until now our discussion of the energy spectra was based 
on states |JV, 0) without bosonic excitations and so far 
the effect of V npp on the spectrum could have even been 
treated without using bosonization. But for the determi- 
nation of the excitation spectrum of H we do need the 
general expression for the matrix elements of V npp between 
the eigenstates of H + V pp as given by (50). For the ac- 
tual calculation we truncate the eigenbasis of Hq + V pp 
for a fixed charge state N c at a certain excitation energy 
and represent H in this shortened basis. After the diag- 
onalization we find to a good approximation the correct 
eigenstates and eigenenergies of H. For the results shown 
in Figs. 7 to 10 we have checked that convergence has been 
reached, i.e., the extention of the considered basis states 
does not lead to a significant change of the spectrum. 

Exemplarily we present the results for the charge state 
N = An. Similar excitation spectra are found for the other 
charge states. In Fig. 7 we show for comparison and in or- 
der to demonstrate the effect of the non forward scattering 
processes the findings for the "standard" theory, i.e., the 



spectrum of H st — Hq +Vf / / as well as the spectrum of 
the full Hamiltonian H — Hq + V pp + V npp for a (6, 6) arm- 
chair nanotube. Thereby a nonvanishing band mismatch 
ea = is assumed. Striking is the partial breaking of 
the huge degeneracies of the "standard" spectrum. Note 
also the lifting of the spin-charge separation when includ- 
ing the non forward scattering processes. To illustrate this 
point we have indicated eigenenergies of H st including c+ 
excitations by arrows in Fig. 7. 

At higher energies a quasi continuum forms in the case 
of the full Hamiltonian H , a feature becoming especially 
apparent for a finite band mismatch. In Fig. 8 the spectra 
of the full Hamiltonian H is shown for ea = 0.3eo. 

As we have already discussed, the importance of non 
forward scattering terms should decrease with increasing 
tube diameter. And indeed the excitation spectrum of the 
full Hamiltonian for a (20, 20) SWNT resembles much 
more the result of the "standard" theory than it is the 
case for a (6,6) SWNT as it can be seen from Fig. 9. 
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Fig. 9. The excitation spectrum for a (6,6) SWNT (black) 
compared to the spectrum of a (20,20) SWNT (grey). The 
effects of the non forward scattering processes are by far less 
pronounced in the latter case. N c — 4m and sa = 0. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 

Number of Eigenvalue x10 4 



situation the meanfield Hamiltonian can essentially be 
recovered by setting all off-diagonal elements in (55) to 
zero. Therefore in [10] the degenerate triplet state can not 
be predicted but twofold degeneracies of the states |T, T), 
||,|) and of |t, 4), ||, T) respectively are found. Moreover 
contrary to our theory in [10] no mixing of the states 
|t|, — ) and | — , tJ.) can occur for sa < J/2, an important 
point regarding the singlet-triplet Kondo effect [19]. 

Moreover also the excitation spectrum shows consid- 
erable differences in both approaches, since the meanfield 
approach misses the formation of the collective electronic 
excitations as the c+ mode, with its dispersion relation 
strongly renormalized by the forward scattering part of 
the Coulomb interaction. 



4.4 Near half-filling 

As we have already seen in Section 3.2, at half-filling 
non-density-density interaction processes become relevant 
which yield tremendously large matrix elements in the 
eigenbasis of Hq + V pp , as a consequence of the function 
A u (x) shown in Fig. 4. Therefore our truncation scheme 
for diagonalizing the total Hamiltonian V npp does not give 
reliable results at half-filling. Investigation of the half- 
filling case is beyond the scope of this work. 



Fig. 10. The excitation spectrum for a (6,6) SWNT obtained 
by diagonalizing the Hamiltonian of the noninteracting system 
Ho (grey) and the full Hamiltonian H — Ho + V pp + Vnpp 
(black). 

It is also interesting to regard the effect of the total 
interaction V pp + V npp on the nanotube spectrum. For this 
purpose, in Fig. 10 the spectrum of Hq describing the 
noninteracting system is compared to the spectrum of the 
full Hamiltonian H, again for a (6, 6) SWNT with vanish- 
ing band mismatch. Of special significance is the strong 
reduction of the number of eigenstates below a certain en- 
ergy if the interaction is "switched on". This can be mainly 
traced back to Vfff which leads to the formation of the 
bosonic c+ excitations with considerably enlarged ener- 
gies. Concerning the transport properties of SWNTs the 
reduction of relevant states plays an important role for the 
occurrence of the power law dependence of various trans- 
port quantities in the case of infinitely long tubes but also 
for the appearance of negative differential conductance in 
highly asymmetric SWNT quantum dots as described in 
[6]. 

4.3 Comparison to the mean field results 

We shortly want to compare the results of the meanfield 
theory by Oreg et al. [10] and our approach. Concern- 
ing the groundstate structure, differences between the two 
works arise for the N c = 4m + 2 charge state. In this 



5 Conclusions 

In summary, we have derived the low energy Hamilto- 
nian for metallic finite size SWNTs including all relevant 
interaction terms, especially the short ranged processes 
whose coupling strength scales inversely proportional to 
the SWNT size. The Hamiltonian of the noninteracting 
system, Ho, together with the density-density part of the 
interaction, V pp , could be diagonalized by bosonization 
and Bogoliubov transformation. Considering only the sit- 
uation away from half-filling, we obtained the spectrum 
of the total SWNT Hamiltonian by exploiting the small 
magnitude of the non-density-density contribution V npp to 
the interaction: we have calculated the matrix elements of 
V npp in a truncated eigenbasis of Hq + V pp and diagonal- 
ized the resulting matrix to obtain the SWNT spectrum 
and the corresponding eigenstates. 

Of special interest, concerning the ground state spec- 
tra, is the formation of a spin 1 triplet for the charge state 
N c = Am + 2, whose existence has clearly been proven in 
the experiments of Moriyama et al. [9]. In the case of a 
band mismatch sa that is small compared to the exchange 
energy J, the spin 1 triplet is the ground state of the sys- 
tem. This finding is interesting since according to a theo- 
rem by Lieb and Mattis [11], only ground states with spin 

or 1/2 are allowed for a ID Hubbard model with next- 
neighbour hopping and no orbital degeneracies. Since our 
SWNT Hamiltonian includes an orbital degree of freedom 
we conclude that scattering processes with respect to this 
degree of freedom are the reason for the finding of a spin 

1 ground state. Additionally we predict for ea < <//2, the 
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mixing of the states |TI,— ) and | — ,TI) with an accom- 
panying energy splitting. The degree of mixing between 
ItJ., — ) and | — , tl) is of importance for the singlet-triplet 
Kondo effect, as discussed in [19]. An experimental confir- 
mation of our findings in the case ea <//2 is still missing, 
but well within reach. 

With regard to the excitation spectrum, the different 
bosonic modes are mixed by the non-density-density in- 
teraction processes V npp . Therefore the spin-charge sep- 
aration is lifted. Moreover we find that the huge degen- 
eracies which are obtained by the "standard" theory that 
retains only forward scattering processes are partially bro- 
ken. This leads to a more and more continuous energy 
spectrum for higher energies. 

A Modelling the interaction potential 

In this Appendix we show how we determine the values of 
the effective ID potentials jj^ ra / mter an( j f ^ ne coupling 

constants u Sr Sf . We start with equation (12) from section 
2.2, 



U™\ ra/inter {x,x')=IS 



a rid v i 



x l P*pF 1 ( r ) ( P*± P F 2 {r')<P± P F 3 (r')<fi P F 4 (r)U(r - r'). (56) 

Using equation (4) in order to reexpress the Bloch waves 
(p P F{r) in terms of p z orbitals, we obtain, 



U, 



intra / inter 
IF] 



(x,x') = 



Nt 



d 2 r ± d 2 r' ± 



x U(r - r') ^2 e~* (Fl ~ F4)flx ~ i(F2 ~ F3)ii * 



R,R' 

x \ x ( r -R-T p )\ 2 \x(r'-R'-T ±p )f. (57) 

Instead of a fourfold sum over the lattice sites only the 
double sum ^2 R R , remains, since the overlap of different 
p z orbitals can be neglected. To proceed we use once more 
that the spatial extention of the p z orbitals is small com- 
pared to all other appearing length scales and therefore re- 
place \\{r — R — r p )| 2 by the delta function 5(r — R—T p ). 
In order to take into account the error induced thereby at 
small distances i ~ i', we replace the Coulomb poten- 
tial by the Ohno potential introduced by equation (8). It 
interpolates between Uo, the interaction energy between 

2 

two p z electrons in the same orbital and A7r e t\r-r'\ ^ or 
\r — r'\ 3> 0. Performing the integration in (57), we ob- 
tain, 



U, 



intra / inter 
IF] 



L 



x e 



L R,R 

i(F 1 -F4)R x -i(F2-F 3 )R' a 



Using (58) together with equation (11), 



U [r][F] (x,x') = - lU$ ra (x,x')(l + nr2r3r4) 
+ U^ix^'^rs+nn) 



(59) 



we arrive at 



uf b =:u+ 



— T 

AN 2 ^ 



i1K (R x -R' x ] 



R.R' 



x[U{R-R') + U{R-R' + T P -T- P )}, (60) 



u W*f =:u f = 

[U(R-R')-U(R -R' + t p - t_ p )] (61) 



L R.R' 



and 



u v«/ =: u f = _L y e -i2K 0{Rx -R' x) 

AN? ^ 



R.R' 



x[U(R-R')-U(R-R' + T P -T- P )]. (62) 

Since in the summations in (60), (61) and (62) only terms 
with R w R' contribute, the number of relevant sum- 
mands scales like the number of lattice sites Nl. Due to 
the prefactor 1/N 2 , u + and u^ b in total scale like 1/Nl- 
Numerical evaluation of the previous three equations leads 
to the values given in table 1. 



B Calculation of the matrix elements 

M [r][F][a] (N,m,N',m',x) 

Using the bosonization identity (39), 

4>r* F (x) = Vra K raF (x)e^^e^^ x \ 

we can separate M^ F ^(N ', m, TV', m! , x) from equation 
(38) into a bosonic and a fermionic part, 

M [r][F][cr](N,m,N',m',x) = 

M [r] [F ] [a] {N,N', x)M {r] [F] [ct] (m,m',x), 

where 

M {l] (N,N',x) = 
(N\ 4 {x)r,l 4 {x)riK h {x)n h K u (x)rj h \N>) (63) 

and 



U(R-R' + t p -t ±p ). (58) M[l] (m,m',x) = (m\ eXWe^iWeXWe^W 



Now we can easily calculate the values of the coupling 
constants u SrSp for the local interactions, given by (20), 



S r Sf 



= l/(2£ 2 ) J J dxdx'U [r]Sr[F]sF {x,x'). 



e i^) e i4>, 3 ^) e ^l(^ e i<t>i^) | m ') . ( 64) 

Improving readability, we have summarized the indices 
rFu by a single index I. 
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B.l The Fermionic part of M [r][F][<T] (N,m,N',m',x) w h ere EjLi^n denotes the sum <j> h + 4> h - <j> h - (f> h and 



First we consider the contribution Mpj (TV, N 1 , x) depend- 
ing on the fermionic configurations N and N'. Using re- 
lation (40) for the Klein factors rj rcr and the definition of 
the phase factor K r(jF (x), equation (41), we obtain 

M m[a] (N,N',x) = 

(2L) 2 ^ JV,W ' +E[r| w Tnn ' M [ ff l ® NN ' M i p ] ( x ) ' 



where E 



r ][ CT ] := e riCT + e r2Cr / - e r3a > - e r4(T . Furthermore 



W MW = ( -l ) ^ 4> - 1 (-')-^ 3> - 1 (-'--- 4 ^3 

x ( _ 1) E^? , - 1 (iV-e. 1CT1 ) 32+2 :<-? ) - 1 (JV) 31 ^ (65) 



Here we use the convention j = + f,+ J,,— 1 = 

1,2,3,4. It turns out that Tjvjv' [r] [c] only depends on the 
scattering types S r and S a . Explicitly with Tjvs r s„ := 

T NN'lr] Sr [<x] S<T > 

T w - =-(-l) 3JV «T+2^i+^T, (66) 

Tjv'6/- = (-l) 3JV «T+ 2Jv k+ w M (67) 

and Tjvs r 5„ = 1 for all other (5 r , Finally the func- 
tion Q jvjv' [r] [F] i x ) yields a phase and is given by 

QNN'[r][F](x) = 

exp [sgn(r 4 F 4 )(N') u + sgn(r 3 F 3 )(;V' - e u ) h 
sgn(r 2 F 2 )(N - e h ) h - sgn(r 1 F 1 )(AT) ;i ] x} . (68) 

B.2 The bosonic part of M[ r ][ F ][<T](N,m,N',m',x) 

The calculation of the bosonic part Mr r irCTr .i(m, m', x) is 
based on expressing the fields i(j> ra - F (x) in equation (64) in 
terms of the bosonic operators djs q , Oj Sq and subsequent 
normal ordering, i.e., commuting all annihilation opera- 
tors djSq to the right side and all creation operators Oj S 
to the left side. In a first step we use the relation 

e i<t>l(x) e i<t>l(x) _ e %4>\{x) e %4>i{x) e [i<t>i(x),i4>l(x)] 

following from the Baker-Hausdorff formula [17], 

e A e B = e A+B e J[A,B] j [A,B]€C, 

to obtain from (64), 
M {l] {m,m',x) = C [t] (x) 



C[l](x) — e ^ l 3^' i ^l 4 ( X ^ e l- i ^2( X )' i ^l 3 ( X )+ i( l>l 4 ( X )] 

x e l~ i 4>l 1 (x)--i4>l 2 (x)+i4>l 3 (x)+i<p f li (x)] 

Applying the Baker-Hausdorff formula once more, we ob- 
tain 

e -iELl(^n(^)+'Ai t n (^)) e 3[ i ^»=l^n( X )' i £n' = l^n'( X )] 

Using the definition of the 0-fields, equation (42), together 
with the transformation between the operators b aq and 
Cbj$ q , equation (33), we get 



i(j) rijF {x) +i4>l aF (x) 



jSq>0 

In terms of A 3 r s a , Bj$ q and Djg q , cf. equations (34), (35) 
and (36), the coefficients X 3 r J F (x) read 



K 5 M = % 



^ e isgn(rF)qx _ e -isgn(rF) q x JJjg^j 



(70) 



By defining 

HtW^-tLi^w (71) 

and again using the Baker-Hausdorff formula, we arrive at 

e -iff n =i{<t>l n {x)+4>] n {x)) = 
e ~12js q> o ^*i] 5q ( x ) a ]s q jS q >o ^ii] q ( x ) a j5q e ~i 52j5 q> a\^ii] q ( x )\ 

such that in total 

•ELl^LWp-'ELl^nW 



m 



rn ) = 



A [q (x) J] FCX>* q (x), m j5q , m' jSq ), (72) 
where we have introduced 



x g 2 E^gX)!^!]"^)! _ ^3) 



m ) is given 



The function F(X,m,m') = (m 
by [6] 

F(X, m, m') = 

(G{m' - m)X m '- m + 0{m - m') (_A*) m " m ') 

(-IAI 2 ) 



i=0 



!(i -)- TTlmax ^miri)^- (j^min ^)! 



(74) 
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where m min / max — min/ max(m, m'). Combining (69) In total this leads to 
and (72) we finally obtain 



M [V[ {m,m!,x) = C[i](x) 



A 



2E >o 7Tr(l-cos(2gx)) 



X A [l] ( X ) II F ft[l] q ( X )> m ^«' m 'jSq)- 



x e 



-2E,>„^(^ + l) S in 2 (^) 



Explicitly, equation (73) yields that Am (a;) only depends 
on the scattering type for the product rF. For S r F /fiwe 
find A[i] s ^ =: As rF = 1 whereas A u is strongly enhanced 



Because of sin 2 (gx) = h (1 — cos(2<7x)) the final result is 



(x) 



iE 9> o^(l-^) S in 2 («x) 



(76) 



leading to an increased importance of non-density-density ml , , . , , , . , „ , , 

°. , , , r £nn . -|— . . i J-he same result is also obtained for all other processes 

interactions at half-fallmg. Due to its relevance we show . ^ S — u 

the detailed calculation of A u in the following. 



B.2.1 Evaluation of A u 



C Regularization of (Nm \ Vf+ b /- 1 Nm) 



As example we calculate ^[ r ] Sr [F| Sf ,[ < T]s < , with (S r , S F , S a ) = As already mentioned in the main text, expression (50) 

(b, f~,f + ), i-e., for [r] = (r, -r, r, -r), [F] = (F, -F, -F, F) for the matrix element (Nm \V Sr s F s a I N'm'} diverges if 

and [a] = (a, a, a, a). It is easily checked that for this v. m jSq - m' jS ^ < 1 and if ^s r S F S„ is N conserving. 

choice S r F = u holds. Before starting with the actual cal- TT , ,, , ■ , , 

, . „ , . , _ . rjSa , \ r Here we show m detail how the matrix element can be 

culation we first determine the coefficients A, „ CT r Ax) for , , ■ , c „. i , T/ 

H[F]M V y properly regularized for the case m — m and Vs r s P s„ = 

the considered case. With equations (70) and (71) we find V f+bf - . We start with equation (50), 



A' 



(x) = 



jSq 

y]biF] f -[o-] f +' 

AJ S I isgn(r n F n )qx r> _ - 



(Nm \ Vf+ bf -\Nm) 



isgn(r n F n )qx 



rFcr J 



dx 



The values for Bjs q , Djs q and A£* are known from the 
Bogoliubov transformation, cf. equations (34) to (36). For 
the different channels jS this leads to 



X [r] b [F] f _ [a] f+ ( X ) 



2^sgn(rF) [T^ 



-c+q 



■ sm(qx), 



\r] b [F] f _[a] f+ W ~ ^ sm{qx), 

(*)=«■ 

Using (73) we get in this case, 

Am (a;) := ^[^h^'^i^^ 1 ^-^ 1 ^] 



x e 



5 [*^I 2 ( x )~ i( t>l 4i ( x )- i <t>l 2 fa)-*^ ( x )] 



x e 



2 



( \\" + q 



(x)| +^=+"(^1 



• (75) 



Improving readability we have again replaced the indices 
rFa by a single index /. With (42) we obtain 



g>0 q 



sgn(rF)q(x^fx) 



V'l 



E 1 

q>Q q 



crr-qi u <Jr-q\ ~ 

isgn(rF)q(x^fx) 



e -2isgn(rF)(N r<T -N r _„)fx 

4 sin 2 (fx) 



x n F (K]f+[F] b [«],- {x),m j5q ,m j5q ). (77) 



In a first step we rewrite the fraction 



-2i SS n(rF)(N rtT -N r _^)f a 

4sin 2 (-£-:r) 



as 



-2iagn(rF)JV ro - ^-x 2isgn (rF ) N r _ a -g- x 

-g^ 2 -i'^ x an dj by using the iden- 



tity 



N 

E 



-iNx 



1 



we transform it into the product of two infinite sums ex- 
tending over the whole Fermi sea, 

-2isgn(r_F)(Ar r „-W r _„)f x 



4sin 2 (fx) 

N rt7 

E 



-22sgn(r.F)r 



E 



■ (78) 



n— — oo 



n = — oc 



An important observation is, that the multiplication with 

e -in q x _ e in qX _ e -»n,x ^ _ g2in,ij recasts the infinite 



-2inx 



into a finite sum, 



N 

E 



TV 



e - 2ma: . (79) 
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We now have a closer look at the coefficients Af 5 , 9 , P1 , , 
which according to (71) are given by 



[jSq 



] s +[F] b [a] s 
1 



(x) 



y[jS ^ e -isgn(rF)qx _ e isgn(rF)qx^j 



+ A* S ^ e is S n ( rF )qx _ e -isgn(rF)qx 

Then because of (79), the product 

n(^> ]bM/ >)) 



jSq 



N r - a 



x e -2^(rF)nix ^ e 2i S gn(r_F)n'f x^ ^ g ^ 



n — — oo 



is a finite sum for *j<5<? — 2- But from (74) we can 

conclude that 

where 0(A 2 ) collects all those terms which contain a factor 

UjSq (^UFU^))^ ^th EjSqhSq > 2. ThuS 



-2isgn(rF)(N ra -N r _ a )fx 



dx- 



4 sin 2 (fa;) 

\jSq 



< 1 II F C>ir] f+ [F] b[ * ]f _ (*). mjtf,, m j4 ,) - 1 I (80) 

is a well defined integral over a finite sum and therefore 
not diverging. On the other hand we find with (78) 



/ 



dx 



e -2iagn(rF)(N ra -N r - a )fx 



4 sin 2 (fa;) 

J dx ^ e -2isgn(rF)n^x ^ ( .2i*nn{iJ 7- ! 



e b K j L y g 

n— — 00 n' — — 00 

AT r<T N r - a mm(N r<T ,N r -„-) 



n— — 00 n'=— 00 



Regularization of the previous expression now is easily 
achieved by subtracting in the previous equation e.g. the 
contribution from below half-filling, such that, 



/ 



e -2isgn(rF)(N rtT -N r - tT )fx 

dx n ; - — ^ = L mm(N r(T , N r _ a ). 



4 sin 2 (fa;) 



(81) 



Combining (80) and (81) we obtain the finite expression, 
(Nm \ V f+b f-\Nm) = u + ^ min(iV rT , N rl ) 



+ JE u+ ^l dx 



rFa 



-2isgn(rF)(N rcr -N r - t ,)%x 



4 sin 2 (fa;) 



n F (5 /+ H.W r {x),m j5q ,m j5q ) - I 



which is equivalent to equation (51) in the main text. The 



m is q ~ m jSq 



— 1 as well 



regularization for the case J2jS q 
as for the matrix elements of the N conserving processes 
Vf-tf and V b f-f+ which are only relevant near half-filling 
can be achieved in a similar way. 
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